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ABSTRACT 

We compare the shapes and intrinsic alignments of galaxies in the MassiveBlack-II 
cosmological hydrodynamic simulation (MBII) to those in a dark matter-only (DMO) 
simulation performed with the same volume (100/i _1 Mpc) 2 3 , cosmological parameters, 
and initial conditions. Understanding the impact of baryonic physics on galaxy shapes 
and alignments and their relation to the dark matter distribution should prove useful 
to map the intrinsic alignments of galaxies from hydrodynamic to dark matter-only 
simulations. We find that dark matter subhalos are typically rounder in MBII, and the 
shapes of stellar matter in low mass galaxies are more misaligned with the shapes of 
the dark matter of the corresponding subhalos in the DMO simulation. At 2 = 0.06, 
the fractional difference in the mean misalignment angle between MBII and DMO 
simulations varies from ~ 28% — 12% in the mass range 10 10,8 — 6.0 x 10 14 /i _ 1 M 0 . 
We study the dark matter halo shapes and alignments as a function of radius, and 
find that while galaxies in MBII are more aligned with the inner parts of their dark 
matter subhalos, there is no radial trend in their alignments with the corresponding 
subhalo in the DMO simulation. This result highlights the importance of baryonic 
physics in determining the alignment of the galaxy with respect to the inner parts of 
the halo. Finally, we compare the ellipticity-direction (ED) correlation for galaxies to 
that for dark matter halos, finding that it is suppressed on all scales by stellar-dark 
matter misalignment. In the projected shape-density correlation (ws+), which includes 
ellipticity weighting, this effect is partially canceled by the higher mean ellipticities of 
the stellar component, but differences of order 30 — 40% remain on scales > 1 Mpc 
over a range of subhalo masses, with scale-dependent effects below 1 Mpc. 

Key words: cosmology: theory - methods: numerical - hydrodynamics - gravita¬ 
tional lensing: weak - galaxies: star formation 


1 INTRODUCTION 

Weak gravitational lensing is a cosmological probe that has 
the potential to address major outstanding cosmological 
problems, such as understanding the connection between 
dark matter and galaxies, the nature of dark energy, and 
exploring possible variations in the theory of gravity on cos¬ 
mological scales (Albrecht et al. 2006; Weinberg et al. 2013). 
Using weak lensing, upcoming surveys such as the Large 
Synoptic Survey Telescope 1 (LSST; LSST Science Collabo- 
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ration et al. 2009), Euclid 2 (Laureijs et al. 2011), and the 
Wide-Field Infrared Survey Telescope 3 (WFIRST; Spergel 
et al. 2015) will to constrain cosmological parameters such 
as the dark energy equation of state to a very high preci¬ 
sion. However, the most basic weak lensing analysis assumes 
that galaxy shapes are randomly aligned, which is not cor¬ 
rect in reality. The galaxy shapes are correlated with each 
other and with the underlying density field. This intrinsic 
alignment of galaxy shapes is an important astrophysical 
systematic that contaminates measurements in weak lens¬ 
ing surveys (Heavens, Refregier & Heymans 2000; Croft & 


2 http://sci.esa.int/euclid/, http://www.euclid-ec.org 

3 http://wfirst.gsfc.nasa.gov 
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Metzler 2000; Catelan, Kamionkowski k Blandford 2001; 
Jing 2002; Hirata k Seljak 2004). For reviews on intrinsic 
alignments and its impact on cosmological parameter con¬ 
straints, see Troxel k Ishak (2014), Joachimi et al. (2015), 
Kirk et al. (2015), and Kiessling et al. (2015). 

Previous studies of intrinsic alignments that made pre¬ 
dictions out to tens of Mpc scales have generally involved 
analytical methods, or iV-body and hydrodynamic simula¬ 
tions of cosmological volume. Analytical modeling involves 
a number of different possible models, with the one receiv¬ 
ing the most attention for elliptical galaxies being the lin¬ 
ear alignment model (Catelan, Kamionkowski k Blandford 
2001; Hirata k Seljak 2004), an extension of it to nonlinear 
scales using the non-linear power spectrum (Bridle k King 
2007; see Blazek, Vlah k Seljak 2015 for a more recent ex¬ 
tension that includes all non-linear terms at the same order), 
and a small-scale extension using the halo model (Schneider 
k Bridle 2010). iV-body simulations have been used to study 
intrinsic alignments by Heymans et al. (2006) where halos 
are populated with galaxies which are stochastically mis¬ 
aligned, Joachimi et al. (2013) with semi-analytic models, 
and others (see references in the review by Kiessling et al. 
2015). However, the physics of galaxy formation is not taken 
into consideration by these analytic models and iV-body sim¬ 
ulations. Recent hydrodynamic simulations of cosmological 
volumes such as the MassiveBlack-II (Khandai et al. 2014), 
the Horizon-AGN (Dubois et al. 2014), and the EAGLE and 
cosmo-OWLS (Velliscig et al. 2015) simulations have made 
it possible to directly study the intrinsic alignments of the 
stellar component of galaxies and their scaling with mass, 
luminosity and distance. 

Using the MassiveBlack-II (MBII) simulation, we pre¬ 
viously studied the shapes of the stellar matter component 
of the galaxies and its alignment with the shape of the dark 
matter component (Tenneti et al. 2014). We found that the 
ellipticities of the stellar components of galaxies compare 
favorably with those in observations, and that the shape of 
the stellar component is misaligned with the dark matter 
component, with average misalignment angles ranging from 
~ 10° to 30°. In Tenneti et al. (2015), we extended this 
study to two-point statistics and explored their dependence 
with mass, luminosity and distance. We found that the in¬ 
trinsic alignments of massive galaxies in the simulation have 
a scaling with transverse distance that is consistent with re¬ 
sults from observational measurements. Large volume hy¬ 
drodynamic simulations are proving to be extremely useful 
in quantifying intrinsic alignment signals. However, they are 
also very computationally expensive. Hence, it would be use¬ 
ful to develop methods to paint the intrinsic alignments of 
galaxies onto iV-body simulations. It is known that the halo 
shapes and orientations in iV-body simulations overestimate 
the intrinsic alignments when compared with observational 
results (Heymans et al. 2006). Hence, a direct comparison 
of intrinsic alignments in hydrodynamic and iV-body simu¬ 
lations is necessary. In this paper, we use the MassiveBlack- 
II dark matter-only (DMO) simulation that has been per¬ 
formed with the same resolution, box size, initial conditions, 
and cosmological parameters as the hydrodynamic simula¬ 
tion (MBII) to compare the intrinsic alignments of galaxies 
in MBII and halos in DMO. 

To predict the intrinsic alignment of galaxies from dark 
matter-only simulations accurately, it is important to have 


an average mapping that statistically determines the ellip- 
ticity and orientation of the stellar matter component of 
the subhalo with respect to that of the dark matter compo¬ 
nent. Here, the radius at which the shape of the dark matter 
subhalo is measured matters, as it has been shown using IV- 
body simulations that the shapes of dark matter component 
change with radius (Allgood et al. 2006; Kuhlen, Diemand 
k Madau 2007; Schneider, Frenk k Cole 2012). Further, it 
is also known that the dark matter halo shapes in hydrody¬ 
namic simulations are different from those in iV-body simu¬ 
lations due to the effects of baryonic physics, which leads to 
rounder shapes (Kazantzidis et al. 2004; Bryan et al. 2013). 
Here, we study the radial dependence of the distribution of 
dark matter halo axis ratios in the MBII and DMO simu¬ 
lations, and compare with the shape of the stellar matter 
component in MBII. We also study the orientation of the 
stellar shape with the shape of dark matter component in 
MBII and DMO measured at different radii. The change 
in the orientation of the dark matter shapes in the subhalos 
measured at different radii has been previously studied using 
iV-body simulations (e.g., Bailin k Steinmetz 2005; Schnei¬ 
der, Frenk k Cole 2012). Using small volume hydrodynamic 
simulations, Bailin et al. (2005), Deason et al. (2011) and 
others (see references in the review by Kiessling et al. 2015) 
found that the inner shape of the dark matter component 
is well aligned with the galaxy orientation. However, it will 
be useful to compare the orientation of stellar shape in the 
hydrodynamic simulation with the shape of the dark matter 
component measured at different radii in the matched sub¬ 
halo of a dark matter-only simulation. This comparison is 
what will enable the mapping of intrinsic alignments in hy¬ 
drodynamic simulations onto dark matter-only simulations. 

This paper is organized as follows. In Section 2, we de¬ 
scribe the details of the simulations used in this study fol¬ 
lowed by the method adopted to determine the shapes of 
galaxies and the definitions of two-point statistics. In Sec¬ 
tion 3, we match the subhalos in both simulations and com¬ 
pare shapes and misalignment angles. A brief discussion of 
the comparison of mass functions and dark matter power 
spectrum in the two simulations is also included here. In 
Section 4, we analyze the radial dependence of shapes and 
orientations of the dark matter component in subhalos. In 
Section 5, we compare the two-point correlation functions 
using shapes defined by the dark matter component in MBII 
and DMO, and the stellar shapes from MBII. Finally, we 
summarize our conclusions in Section 6. 

2 METHODS 
2.1 Simulations 

In this study, we use the MassiveBlack-II hydrodynamic 
(MBII) and Dark Matter Only (DMO) Khandai et al. (2014) 
simulations. MassiveBlack-II is a state-of-the-art high res¬ 
olution, large volume, cosmological simulation of structure 
formation. These simulations have been performed in a cubic 
periodic box of size 100/i _1 Mpc on a side using P-GADGET, 
which is a hybrid version of the parallel code, GADGET2 
(Springel, Di Matteo k Hernquist 2005) that has been up¬ 
graded to run on Petaflop-scale supercomputers. The total 
number of dark matter particles in both the simulations is 
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Figure 1 . Snapshot of the MBII (left) and DMO (right) simulations in a slice of thickness 2h -1 Mpc at z = 0.06. The shading represents 
the density distribution of dark matter, and the red lines show the alignment of stellar shapes in MBII. The length of the lines is 
proportional to the size of major axis of the ellipse representing stellar shape. 



Figure 2. Snapshot of a massive halo (~ 10 14 h —1 M©) in the MBII (left) and DMO (right) simulations, showing the density distribution 
of dark matter at z = 0.06. The blue and green circles show the virial radii of the subhalos centered at the location with highest density. 
In both the panels, ’x’ and indicate the locations of the halo centers in the MBII and DMO simulations, respectively. 


1792 3 with an equal initial number of gas particles in the 
hydrodynamic run. The cosmological parameters are chosen 
according to WMAP7 (Komatsu et al. 2011), with amplitude 
of matter fluctuations as = 0.816, spectral index n s — 0.96, 
matter density parameter = 0.275, cosmological con¬ 
stant density parameter Qa = 0.725, baryon density param¬ 
eter Qb = 0.046 (in MBII), and Hubble parameter h — 0.702. 
Table 1 shows the box size (Lb ox ), force softening length 
(e), total number of particles including dark matter and gas 


(Apart), mass of dark matter particle (tudm) and initial mass 
of gas particle (m gas ) for the two simulations. MBII includes 
the physics of a multiphase interstellar medium (ISM) model 
with star formation (Springel &; Hernquist 2003), black hole 
accretion and feedback (Springel, Di Matteo & Hernquist 
2005; Di Matteo et al. 2012) in addition to gravity and 
smoothed-particle hydrodynamics (SPH). Radiative cooling 
and heating processes are included (as in Katz, Weinberg & 
Hernquist 1996), as is photoheating due to an imposed ioniz- 
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Table 1 . Simulation parameters: box size (Lb ox )> force soften¬ 
ing length (e), number of particles (A/part), mass of dark matter 
particle (ttidm) and mass of gas particle (ra gas ) 


Parameters 

Hydrodynamic 

(MBII) 

Dark Matter Only 
(DMO) 

L b ox (A -1 Mpc) 

100 

100 

e (A -1 kpc) 

1.85 

1.85 

X part 

2 x 1792 3 

1792 3 

m DM 

1.1 x 10 7 

1.32 x 10 7 

m g as(A _1 M 0 ) 

2.2 x 10 6 

0 


ing UV background. Further details about the MBII simula¬ 
tion can be found in Khandai et al. (2014). The dark matter- 
only simulation, DMO, is performed with the same volume, 
resolution, cosmological parameters and initial conditions 
as in MBII. Figure 1 shows snapshots of the dark matter 
distribution in a slice of 2/i _1 Mpc thickness at z = 0.06 in 
both the MBII and DMO simulations. As expected, the dark 
matter distributions appear to be nearly identical on large 
scales. For further comparison, we also show the dark mat¬ 
ter distribution in an isolated halo in Figure 2. Here, we can 
clearly see small differences in the distribution of dark mat¬ 
ter within the dark matter halo virial radius which, as we 
shall see, can lead to changes in the shapes and orientations 
of the halos and subhalos. 

The halo catalogs in the simulations are generated us¬ 
ing the friends of friends (FoF) halo finder algorithm (Davis 
et al. 1985) with a linking length of 0.2 times the mean inter¬ 
particle separation. The subhalo catalogs are generated us¬ 
ing the SUBFIND code (Springel et al. 2001) on the halo cata¬ 
logs. Here subhalos are locally overdense, self-bound groups 
of particles within the halo. Groups of particles are iden¬ 
tified as subhalos if they have at least 20 gravitationally 
bound particles. When performing our analysis, however, we 
use only subhalos from the hydrodynamic simulation with 
at least 1000 dark matter and star particles based on the 
convergence test in Tenneti et al. (2014). 


2.2 Shapes of dark matter subhalos 

In this section, we describe the measurements of dark mat¬ 
ter and stellar matter component shapes in subhalos. We 
model these shapes as ellipsoids in three dimensions using 
the eigenvalues and eigenvectors of the reduced inertia ten¬ 
sor. The projected shapes are calculated by projecting the 
halos and subhalos onto the XY plane and modeling the 
shapes as ellipses. In 3D, the eigenvectors of the inertia ten¬ 
sor are e a ,e b ,e c with corresponding eigenvalues X a ,X b ,X c , 
where A a > A& > A c . The eigenvectors represent the prin¬ 
cipal axes of the ellipsoid, with the lengths of the principal 
axes (a, b , c) given by the square roots of the eigenvalues 
(a/A a, y/Xc). The 3D axis ratios are defined as 


In 2D, the eigenvectors are e a ,e b with corresponding 
eigenvalues A' a , X' b , where X' a > X b . The lengths of the major 
and minor axes are a' = y/X^ and b' = y/X^ with axis ratio 
q 2 d = b'/a'. 

The shapes are determined using an iterative method 


based on the reduced inertia tensor: 


Iij — 


E„ m n 


where 


En ™„ 

= 


( 2 ) 

( 3 ) 


This definition of inertia tensor gives more weight to par¬ 
ticles that are closer to the center which is desirable since 
it eliminates the bias due to loosely bound particles present 
in the outer regions of the subhalo. In the iterative method, 
we first determine the axis ratios by the standard defini¬ 
tion of the reduced inertia tensor using all the particles of 
a given type in the subhalo. Keeping the enclosed volume 
constant (as in Schneider, Frenk X Cole 2012), we rescale 
the lengths of the principal axes of ellipsoids accordingly. 
After this rescaling, we determine the shapes again, discard¬ 
ing particles outside the ellipsoidal volume. This process is 
repeated until convergence is reached. Our convergence cri¬ 
terion is that the fractional change in axis ratios must be 
below 1%. We considered the impact of using different def¬ 
initions of the inertia tensor on the distributions of shapes 
and the intrinsic alignment two-point correlation functions 
in a previous study (Tenneti et al. 2015). 


2.3 Misalignment angle 

To study the relative orientation between the shapes de¬ 
fined by the dark matter and stellar matter components in 
subhalos, we compute the probability distributions of mis¬ 
alignment angles as in Tenneti et al. (2014). If ea a and e ga 
are the major axes of the shapes defined by the dark mat¬ 
ter and stellar matter components, respectively, then we can 
define the misalignment angle by 

6m = arccos(|e da • e ga \). (4) 

In previous work (Tenneti et al. 2014, 2015), we studied 
the distribution of misalignment angles between the shapes 
of the stellar component and dark matter component in 
the subhalos of MBII. Here, our aim is to study the mis¬ 
alignment angles between the shapes of stellar component 
in MBII and the dark matter component in the correspond¬ 
ing subhalo of the DMO simulation. Since the simulations 
have been performed with the same initial conditions, this 
is a well-defined comparison, and the resulting distributions 
will be helpful to paint galaxies onto iV-body simulations by 
providing the necessary probability distributions from which 
to draw the galaxy orientation. 


2.4 Two-point statistics 

To quantify the intrinsic alignments of galaxies with the 
large-scale density field, we use the ellipticity-direction (ED) 
and the projected shape-density (ic^+) correlation functions. 
The ED correlation quantifies the position angle alignments 
of galaxies in 3D, while the projected shape correlation func¬ 
tion can be used to compare against observational measure¬ 
ments that include the 2D shape of the galaxy. 

The ED correlation function cross-correlates the orien¬ 
tation of the major axes of subhalos with the large-scale 
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density field. For a subhalo centered at position x with ma¬ 
jor axis direction e a , let the unit vector in the direction of a 
tracer of the matter density field at a distance r be r = v/r. 
Following the notation of Lee et al. (2008), the ED cross¬ 
correlation function is given by 

u(r) = (| e a (x) -r(x) | 2 } - 1 (5) 

which is zero for galaxies randomly oriented according to a 
uniform distribution. 

The matter density field can be represented using either 
the positions of dark matter particles (in which case the 
correlation function is denoted by the symbol c us) or the 
positions of subhalos (in which case it includes a factor of 
the subhalo bias, and is simply denoted a;). Here, we only 
use uj§ to eliminate the effect of subhalo bias. 

The projected shape correlation functions are computed 
to directly compare our results from simulations with obser¬ 
vations. Here, we follow the notation of Mandelbaum et al. 
(2006) to define the galaxy-intrinsic shear correlation func¬ 
tion, £s+(r p ,n) and the corresponding projected two-point 
statistic, u> 5 +. Here, r p is the comoving transverse separa¬ 
tion of a pair of galaxies in the XY plane and n is their 
separation along the Z direction. 

If 0 . 2 d — b'/a' is the axis ratio of the projected shape 
of the dark matter or stellar matter component of a sub¬ 
halo, and 0 is the major axis position angle with respect 
to the reference direction (tracer of the density field), the 
components of the ellipticity are given by 

(e+, e x ) = 7 —%^ [cos ( 20 ), sin ( 20 )], ( 6 ) 

1 T #2 d 

where e+ refers to the radial component and e x is the com¬ 
ponent at 45°. The matter-intrinsic shear correlation func¬ 
tion cross-correlates the ellipticity of galaxies with the mat¬ 
ter density field, 

p , ^ S+D - S+R 

fe+(r p ,n)= + RR + (7) 

where s+ represents the “shape sample” which is selected 
on the basis of a threshold or binning in subhalo mass. The 
dark matter particles used to trace the density field form a 
“density sample” denoted by D. S+D is the following sum 
over all shape sample vs. dark matter particle pairs with 
separations r p and n: 

S+D= £ s±+ 2 , ( 8 ) 

where e+(j\i) is the + component of the ellipticity of a 
galaxy (j) from the shear sample relative to the direction of 
a tracer of density field (i) selected from the density sample. 
Here, 7Z = (1 — e 2 ms ) is the shear responsivity that converts 
from distortion to shear (Bernstein X Jarvis 2002), with e rm s 
being the RMS ellipticity per component of the shape sam¬ 
ple. S+R is defined by a similar equation for the correlation 
of the data sample with a random density field distribution 
to remove observational systematics in the shear estimates, 
so we neglect this term here. The RR term in Eq. 7 is given 
by the expected number of randomly-distributed points in a 
particular (r p ,n) bin around galaxies in the shape sample. 

The projected shape correlation function, (r p ) is 



Figure 3. Comparison of the halo mass function in MBII and 
DMO simulations at z = 1.0 and z = 0.06. 

now given by 

/ +n max 

| 5 +(r p ,n)dn (9) 

-n max 

We calculated the correlation functions over the whole 
length of the box (100/i - 1 Mpc) with n max = 50/i -1 Mpc. 
The projected correlation functions are obtained via direct 
summation. 

3 RESULTS 
3.1 Mass function 

In Figure 3, we compare the mass function of halos in MBII 
and DMO simulations at redshifts z = 1 and z = 0.06. At 
both redshifts, there are differences in the dark matter halo 
mass function in MBII. There are ~ 10 — 20% more halos in 
the DMO simulation at low masses (10 9 —10 12 /i _ 1 M©). This 
fraction decreases at higher masses. Velliscig et al. (2014) 
found similar results using the OWLS simulation with re¬ 
solved halos in the mass range 10 12 — 10 15 (based 
on o an )- However, Velliscig et al. (2014) find an even 
larger suppression in their halo mass function, possibly due 
to a stronger AGN feedback model. This suppression of mass 
function can be explained by the decrease in the FOF mass 
of halos in MBII. To illustrate further, we plotted the median 
of the fractional difference in the FOF mass of the matched 
halos of MBII and DMO simulation in Figure 4. We find that 
the halo mass is smaller in the hydrodynamic simulation by 
~ 10 — 20 % in the low mass range. 

3.2 Matter correlation function and power 
spectrum 

Figure 5 shows a comparison of the 3D dark matter two- 
point correlation function and power spectrum at z = 1.0 


© 0000 RAS, MNRAS 000 , 000-000 











6 


Tenneti et al. 



Figure 4. Median and scatter (defined using the 16 th and 
84 th percentiles) of the fractional difference in the FOF mass of 
matched halos of MBII and DMO simulations at z = 0.06. 


and z — 0.06 in the MBII and DMO simulations. To cal¬ 
culate the two-point correlation function, we randomly sub¬ 
sampled 4 x 10 5 dark matter particles and tested that the 
correlation function has already converged using a smaller 
subset. The power spectrum is calculated by taking the 
Fourier transform of the two-point correlation function, 

P(k) = 4 tt / £(r)r— ^ —-dr. (10) 

Jo kr 

From the left panel of the figure, we observe that £(r) is 
larger at small scales (< l/i _1 Mpc) in the MBII simulation, 
presumably due to baryonic physics. At these scales, the 
correlation function is larger in MBII by a factor of ~ 10%, 
and the discrepancy is larger at lower redshift. At interme¬ 
diate scales (above ~ l/i _1 Mpc), we find that the ratio of 
the correlation functions approaches unity with no signifi¬ 
cant redshift dependence. The enhancement of the correla¬ 
tion function at small scales is compensated by a decrease 
in clustering at scales comparable to the box size. 

The right panel of Fig. 5 shows a comparison of the 
matter power spectra, demonstrating that the dark matter 
power spectrum in the hydrodynamic run is enhanced by 
~ 10% at k > 10/i/Mpc, with the effect again being stronger 
at low redshift. In contrast to our results, van Daalen et al. 
(2011) found that the dark matter power spectrum is sup¬ 
pressed at intermediate scales in the OWLS simulation that 
included AGN feedback. The purple band in this plot in¬ 
dicates the target precision of ~ 1 per cent accuracy in 
predictions of the matter power spectrum on scales of 0.1 h 
Mpc -1 < k < 10 h Mpc -1 (Huterer & Takada 2005; Hearin, 
Zentner & Ma 2012). This level of accuracy in the theoret¬ 
ical predictions is necessary to avoid systematic errors in 
constraining cosmological parameters. Zentner et al. (2013); 
Eifler et al. (2014) and Mohammed & Seljak (2014) dis¬ 
cuss approaches to mitigate the potentially large differences 


in the matter two-point correlations and reduce the uncer¬ 
tainty to the necessary level. 

3.3 Matching subhalos in MBII and DMO 

To make a fair comparison between intrinsic alignments of 
galaxies and halos in the MBII and DMO simulations, we 
matched the subhalos in both simulations. The subhalos are 
matched using the unique ID’s of the dark matter particles 
in both the DMO and MBII simulations. For a given subhalo 
in the MBII simulation with at least 1000 dark matter and 
star particles, we identify the subhalo in the DMO simula¬ 
tion that has the highest number of dark matter particles 
with the same IDs as in the MBII subhalo. This subhalo 
in the DMO simulation is linked to the subhalo in MBII if 
the fraction of common dark matter particles is greater than 
50% of the total number of dark matter particles in the sub¬ 
halos of each simulation. Figure 6 shows the matched frac¬ 
tion of subhalos as a function of mass. We observe that the 
fraction of matched subhalos increases with mass and ap¬ 
proaches unity. For the rest of this paper, we will only con¬ 
sider subhalos with mass greater than 10 lo ‘ 8 /i _1 M 0 , where 
the matched fraction exceeds 90%. 

The decrease in the fraction of matched subhalos at 
lower masses is due to an increase in the number of satellite 
subhalos. The fraction of central subhalos matched even in 
the lower mass bin, 10 lo ' 8_11 '°/i _1 M©, is ~ 99%. As illus¬ 
trated in Figure 2, there are small but visible differences in 
the smallest subhalo population in a matched halo of mass 
~ 10 14 h~ 1 Mo (for example, some missing or merged sub¬ 
halos in MBII). We also attempted to match subhalos using 
only the 50 innermost dark matter particles, and obtained 
results that were largely consistent with the method dis¬ 
cussed earlier. However, for a small fraction of subhalos, this 
new method would cause us to link subhalos in the two sim¬ 
ulations with very different masses. For example, consider 
the distribution of dark matter in the MBII and DMO sim¬ 
ulations in Figure 2. The blue and green circles show the 
virial radius of the central subhalo and a nearby satellite 
subhalo respectively. We can clearly see that the center of 
the satellite subhalo in MBII is closer to the center of the 
central subhalo in DMO simulation, which leads to a false 
match using the 50 closest particles. Hence, we use all the 
dark matter particles in the subhalo to obtain a consistent 
match. For similar reasons (such as changes in the location 
of density peaks), it is not possible to consistently match the 
satellite subhalos in low mass halos, which is the motivation 
for our adoption of a lower mass limit of 10 10 ' 8 /i _1 Mq. 


4 SHAPES AND MISALIGNMENT ANGLES 

In this section, we investigate the change in the shape and 
orientation of the subhalos with the distance from their cen¬ 
ter, which is defined as the location of the most bound par¬ 
ticle in the subhalo. This measurement is necessary to un¬ 
derstand whether the shape of the stellar matter component 
traces the (inner) shape of the dark matter component. We 
show the radial dependence of the shape of the dark matter 
component in both the MBII and DMO simulations, and 
compare the orientations at various radii against that of the 
stellar component in MBII. 
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Figure 5. Fractional difference between the dark matter two-point correlation functions (left) and power spectra (right) in DMO and 
MBII simulations at z = 1.0 and z = 0.06. The shaded regions represent a deviation within ±1%. 



Figure 6. Fraction of subhalos matched at z = 1.0, 0.3, and 0.06 
as a function of the total subhalo mass in the MBII simulation. 

4.1 Radial dependence of shapes 

To calculate the shape of the dark matter component in the 
subhalo within a given radius, we start with all dark matter 
particles inside the spherical volume at a given distance from 
the center and compute the shape using the iterative reduced 
inertia tensor (Eq. 2). By using the reduced inertia tensor 
for shape calculation, we still provide more weight to the 
particles in the inner region of the subhalo, so our resulting 
shapes should be considered as shapes within a radius rather 
than at that radius. Hence, the effect of baryonic physics on 
the dark matter halo shapes in MBII will be evident even 


for shapes measured out to large radii. For the analysis of 
shapes, we only consider subhalos with at least 500 particles 
remaining after the last iteration, excluding fewer than 1% 
of galaxies. We will compare the distribution of shapes and 
alignments of the dark matter component calculated within 
a radial distance of 0.2R 2 oo, O.6R200 and l.Oifeoo- Here U 200 
is the radius within which the spherically-calculated average 
density of matter is equal to 200 times the critical density. 

The radial dependence of dark matter halo shapes has 
previously been studied using iV-body simulations (e.g., All¬ 
good et al. 2006; Schneider, Frenk & Cole 2012) and small- 
volume hydrodynamic simulations (e.g., Kazantzidis et al. 
2004; Abadi et al. 2010; Deason et al. 2011). As in our work, 
the method used to measure the shapes in Allgood et al. 
(2006), Deason et al. (2011), and Schneider, Frenk &; Cole 
(2012) is based on the iterative reduced inertia tensor. How¬ 
ever, (Kazantzidis et al. 2004) measure shapes using only 
particles within different radial bins with the reduced inertia 
tensor, and Abadi et al. (2010) fit ellipsoids to the positions 
of dark matter particles along the isopotential contours at 
different radii. Hence, the comparison with some of these 
previous studies is only qualitative. 

In Figure 7, we show the normalized histograms of the 
dark matter subhalo axis ratios in MBII and DMO calcu¬ 
lated with different maximum radii, as well as the axis ratios 
of the total stellar matter component of MBII. Through¬ 
out this section, we use three mass bins: 10 lo ' 8_11 ' 5 /i _1 M Q , 
IQ 11 5—13 °/i _ 1 A4© and > lO^/T^M© (Ml, M2, and M3). 
Comparing the dark matter subhalo axis ratio distributions 
in MBII and DMO simulations, we observe that for a given 
maximum radius, the shapes measured in the hydrodynamic 
run are more spherical. This finding is in agreement with 
results from previous studies using smaller-volume hydro- 
dynamic simulations (Kazantzidis et al. 2004; Abadi et al. 
2010; Deason et al. 2011). We also find that the axis ratios 
increase as we go to larger radii, which means that the shape 
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Figure 7. Normalized histograms of 3D axis ratios of the dark matter component in matched subhalos in the DMO and MBII simulations, 
with the shapes measured at different radii (O. 2 R 200 , 0. 6 R 200 and 1. OR 200 ) and also the stellar matter component in MBII. The columns 
indicate different mass bins, while the top and bottom rows are for q (b/a) and s (c/a), respectively. 



Figure 8. Median of the dark matter subhalo axis ratios (left: q , right: s) in the MBII (dashed lines) and DMO (solid lines) simulations 
in different mass bins, plotted against the distance to which the shape is measured, at z = 0.06. 
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Figure 9. Median and scatter (defined using the 16 th and 84 th percentiles) in the distribution of the axis ratios (3D and 2D) of the 
stellar matter component in MBII plotted against the axis ratio of the shape of matched subhalo in the DMO simulation, with the shape 
measured within different radii (O. 2 R 200 , O. 6 R 200 and 1.0i?20o)- Left: q (3D); middle: s (3D); right: q 2 d ( 2 D). 




10 s 



0.2 0.4 0.6 0.8 1.0 
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Figure 10. Contour plots of 3D axis ratios of the dark matter component in matched subhalos in the DMO (left) and MBII (middle) 
simulations, and the stellar matter component in MBII (right), in the mass bin M2 at z = 0.06. 


of the dark matter component is more flattened in the in¬ 
ner regions of the subhalo, again in agreement with previous 
findings (Deason et al. 2011; Schneider, Frenk & Cole 2012). 

To illustrate this effect further, we plot the median axis 
ratios against the radius within which the shape is mea¬ 
sured in Figure 8, in more narrowly-defined mass bins. From 
the plot, we can see that for a given mass bin, the median 
axis ratios are higher in MBII. The median values of q and 
s can be compared against those from Schneider, Frenk & 
Cole (2012), and qualitatively we reproduce their trend that 
they decrease at small distances from the subhalo center 
and at increasing mass. However, this increase in axis ratios 
is smaller in MB-II than in the DMO simulation. Also, at 
higher masses, the increase in the median axis ratio with 
radius is milder than at lower masses. Using smaller-volume 
hydrodynamical simulations, Abadi et al. (2010) found that 
the halo axis ratios are independent of radius, and Deason 
et al. (2011) found that the shapes of dark matter halos are 
slightly more oblate in the inner regions. These differences 
are most likely due to the absence of stellar and AGN feed¬ 
back in those studies, unlike in MBII. 

From Figure 7, we see that the axis ratio histograms 
for the dark matter subhalos in the DMO simulation do not 
follow those for the galaxy stellar components in MBII, at 
any radius. This result can be seen more clearly in Figure 9, 
where we plot the median and scatter in the stellar matter 


axis ratio distributions in MBII as a function of the DMO 
dark matter q or s value, for different radii. We see that the 
scatter in the distribution of axis ratios is large for all radii, 
so there is no advantage in using the inner shape of the dark 
matter subhalo in dark matter-only simulations to predict 
the shape of the stellar component in MBII. For this reason, 
we only consider the dark matter subhalo axis ratios using 
ah particles in our analysis of two-point statistics. 

In Figure 10, we show the contour plot of q versus s 
for the dark matter shape in DMO and MBII and stellar 
shape in MBII in the mass bin M2 (lO 11 ' 5-13 ' 0 /! -1 ./^©). The 
contour plots indicate that the galaxy shapes are more pro¬ 
late compared to the shapes of the dark matter component 
in MBII. We can use the triaxiality parameter, T as 1 _ q 2 
(Franx, Illingworth & de Zeeuw 1991), to quantify the pro- 
lateness or oblateness of the shape. Large (small) values of 
T imply that the shape is more prolate (oblate). In the mass 
bin M2, the mean triaxiality of the stellar component shapes 
is 0.562 zb 0.003, while for dark matter shapes in MBII, 
the mean triaxiality is 0.538 ± 0.002. The triaxialities are 
larger for the dark matter shapes in the DMO simulation 
with a mean value of 0.600 db 0.002, meaning that the dark 
matter shapes are more oblate (prolate) in hydrodynamic 
(dark matter-only) simulations. This conclusion is consis¬ 
tent with results from previous comparisons performed us¬ 
ing simulations of smaller volume (Berentzen &; Shlosman 
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2006; Gustafsson, Fairbairn & Sommer-Larsen 2006; Abadi 
et al. 2010; Deason et al. 2011). 

4.2 Misalignment angles 

In Fig. 11, we show the normalized histograms of the mis¬ 
alignment angles (Eq. 4) between the 3D shape defined by 
the stellar component of subhalos in MBII and the shape de¬ 
fined by the dark matter component in MBII and DMO (for 
mass bins Ml, M2 and M3). From the plot, we can see that 
for M3, the alignment of stellar matter in MBII with the 
dark matter component in DMO (purple curve) is similar 
to the alignment with the dark matter component in MBII 
(dashed black curve). This result implies that the shapes of 
dark matter components in the matched subhalos of MBII 
and DMO have similar orientations at high mass, which is 
also shown directly by the green line. Fig. 12 shows a similar 
result using the projected shapes, with slightly smaller mis¬ 
alignment angles. However, for Ml and M2, it is clear that 
the stellar matter in MBII is better aligned with the MBII 
dark matter subhalo than with the corresponding subhalo 
in DMO. 

To check for a connection between galaxy shapes and 
misalignment angles, we consider the triaxiality parameter, 
T. In Figure 13, we plot the mean misalignment angles be¬ 
tween the shapes of stellar matter in MBII with the dark 
matter shapes in DMO and MBII simulations as a func¬ 
tion of triaxiality. From the figure, we can see that as T 
increases for the stellar shapes and dark matter shapes in 
MBII, the mean misalignment angles decrease. This means 
that for stellar and dark matter shapes in MBII which are 
more prolate, the alignment between the shapes of stellar 
and the dark matter components is stronger than for those 
with more oblate shapes. Similarly, for the more prolate dark 
matter shapes in the DMO simulation, the alignment with 
the stellar component in MBII is closer. However, the mean 
misalignment angles decreases by only ~ 27% while going 
from T < 0.33 to T > 0.66 in DMO, which is less than 
the decrease of ~ 45% for stellar shapes in MBII. When us¬ 
ing galaxy shapes, it is tempting to interpret more oblate 
shapes as relating to disk galaxies, and thus inferring that 
disk galaxy shapes are more misaligned with their host dark 
matter halos than elliptical galaxies. However, we defer a 
more detailed comparison of the morphologies of galaxies 
and the connection to misalignment angle distributions in 
future work. 


4.3 Radial dependence of misalignment angles 

Next, we investigate radial trends in the orientation of the 
shape of the stellar component in MBII with respect to 
the dark matter component at different radii in MBII and 
DMO. Histograms of the misalignment angles when defining 
the dark matter halo shape within various radii are shown 
in Figure 14. For comparison, we also show the prediction 
for a purely random distribution of misalignment angles in 
3D. From the plots, we see that in the MBII simulation, 
the alignment of the stellar component with the dark mat¬ 
ter subhalo shape increases as we limit ourselves to smaller 
radii within the dark matter subhalo. This is expected, as 
the stellar matter is coupled to the distribution and shape 


of the dark matter in the inner regions of the subhalo. How¬ 
ever, there is no corresponding trend in the DMO simulation, 
where the misalignment angle distributions have very little 
dependence on the maximum radius. 

From this plot, we can conclude that it is not advanta¬ 
geous to use the inner shape of the dark matter subhalo in 
a dark matter-only simulation when trying to define mock 
galaxy shapes and alignments. Since the distribution of mis¬ 
alignment angles with respect to the total dark matter sub¬ 
halo orientation is not consistent with a uniform random 
distribution, we can still use these distributions in differ¬ 
ent mass bins to assign shapes and orientations to galaxies 
painted onto IV-body simulations. 


5 INTRINSIC ALIGNMENT TWO-POINT 
CORRELATION FUNCTIONS 

In this section, we compare the intrinsic alignments in MBII 
and DMO simulations by analysing the two-point statis¬ 
tics, ellipticity-direction (ED) and the projected density- 
shape (w$+) correlation functions defined in Section 2.4. In 
Figure 15, we compare the ED correlation for the shapes 
of the dark matter subhalo and stellar matter component 
in MBII with that for the dark matter subhalos in the 
DMO simulation, mass bins Ml (10 lo ' 8_11 ' 5 /i _1 M©), M2 
( 10 ii-s-is-o^-i^©), and M3 (> 10 13 0 /i- 1 M©). The ED 

correlation function for the dark matter subhalos in MBII is 
comparable to that in the DMO simulation on large scales. 
However, on small scales, we observe a tens of percent de¬ 
crease in the ED correlation function for the MBII simula¬ 
tion. This finding is possibly due to the dark matter subhalos 
in MBII being rounder than in DMO, which implies that the 
dark matter particles assume a distribution within the ha¬ 
los that is closer to spherical, reducing the ED correlation 
function at smaller scales. 

The ED correlation computed using the shape of the 
stellar matter component in MBII is smaller than that of 
the dark matter component in DMO (by tens of percent, 
on all scales). This result is due to the misalignment of the 
stellar shape with the dark matter subhalo shapes that are 
determined by the density field. This ratio is relatively scale- 
independent in the higher mass bins M2 and M3, unlike the 
ratio of ED for the dark matter subhalos in MBII vs. in 
DMO. On 0.5 — 5/i _1 Mpc scales, the fractional difference of 
this ratio is on average ~ 72%, 53% and 29% in the mass 
bins Ml, M2, and M3, respectively. 

In Fig. 16, we compare the projected shape-density cor¬ 
relation function, w§+, for the dark matter and stellar mat¬ 
ter shapes in MBII with that of the dark matter shape in 
DMO. Here, ws+ for the dark matter shape in DMO is higher 
than the other correlations at all values of r p and for all mass 
bins. This finding is expected when comparing the DMO and 
MBII dark matter subhalos, since the alignment of the dark 
matter shape with the density field is similar in both simu¬ 
lations, but the dark matter subhalos in MBII are rounder, 
reducing ws+ for MBII. If we consider the ws+ computed us¬ 
ing the shapes of the stellar matter components, we observe 
that at small scales, it is close to the w§+ of the dark mat¬ 
ter shapes in the DMO simulation. This similarity derives 
from the compensation of two competing effects: the stellar 
shapes are more misaligned with the density field, lowering 
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Figure 11. Misalignment angle distributions for the 3D shapes of the dark matter component in matched subhalos of DMO and MBII 
simulations with the stellar matter component in MBII, in mass bins Ml, M2 and M3 at z = 0.06. Also shown (green line) is the 
histogram of misalignment angles between the shapes of dark matter subhalos in MBII and DMO. 





Figure 12. Misalignment angle distributions for 2D shapes of the dark matter component in matched subhalos of DMO and MBII 
with the stellar matter component in MBII, in mass bins Ml, M2 and M3 at z = 0.06. Also shown (green line) is the histogram of 
misalignment angles between the shapes of dark matter subhalos in MBII and DMO. This figure is simply the 2D version of Fig. 11. 


w $+, but the stellar shapes have a larger ellipticity, raising 
w $+. However, the correlation function computed using stel¬ 
lar component shapes is still ~ 30 — 40% smaller than that 
using DMO subhalos for 0.5 < r p < 5/i _1 Mpc. Due to the 
limited size of the simulation volume, the uncertainties are 
large beyond ~ 5/i _1 Mpc. At intermediate scales, we can 
also see the transition from the 1-halo term to the 2-halo 
term for the stellar component, whereas for the dark matter 
component, it is not clearly evident. To illustrate further, 
we also show the ratio for the 1-halo terms corresponding to 
w§+ of the shapes of dark matter and stellar component in 
MBII. In the lowest mass bin, which has a fairly equal mix 
of central and satellite subhalos, we observe a change in the 
shape of the 1-halo term at ~ l/i _1 Mpc for the stellar com¬ 
ponent. As the stellar component is more aligned with the 
inner regions of the dark matter distribution in the subhalo, 
the 1-halo term drops more significantly in the intermediate 
scales in comparison to that of the dark matter component. 

For a more direct understanding of the effects of align¬ 
ment versus different shape distributions, we compute a 
new statistic, ws+, which is defined as the projected shape- 
density correlation function without ellipticity weighting 


(i.e., setting |e| = 1 for all objects in the shape sample). The 
results are shown in Figure 17. For the stellar component in 
MBII, w$+ is smaller than that of the dark matter compo¬ 
nent in DMO at all scales. For 0.5 < r p < 5/i _1 Mpc, the 
fractional differences in the correlation function are ~ 60%, 
49% and 38% in the mass bins Ml, M2, and M3, respec¬ 
tively. Clearly, since the shapes have been normalized to the 
same value, this must be due only to misalignments between 
galaxy shapes and the density field. Similarly, the correla¬ 
tion functions for the dark matter subhalo shapes agree on 
large scales, but that for MBII is smaller than that for DMO 
on small scales, due to the rounder shapes of the dark matter 
subhalos (just as was seen for the ED correlation). 


6 CONCLUSIONS 

In this paper, we carried out a comparison of halo and sub¬ 
halo populations in the MassiveBlack-II (MBII) simulation, 
and a corresponding dark matter-only simulation (DMO) 
with the same cosmology, resolution, volume, and initial 
conditions. First, considering basic halo properties, we com- 
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Figure 13. Mean of the 3D misalignment angles between the 
shape of the dark matter component in MBII and DMO with the 
shape of stellar component in MBII as a function of the triaxiality 
parameter, T (of the stellar matter component in MBII and the 
dark matter component in MBII and DMO), in the mass range, 
10 10 -8 _ g q x 10 14 H~ 1 Mq at z — 0.06. The purple line shows 
the mean misalignment angle between the shape of dark matter 
component in MBII with the stellar component in MBII plotted 
against the triaxiality of the shape of stellar component in MBII. 
Similarly, the green and black lines show the mean misalignment 
angles between the shapes of dark matter component in MBII 
and DMO with the stellar component in MBII plotted against 
the triaxialities of the shapes of dark matter component in MBII 
and DMO respectively. 


pared the halo mass function of the MBII simulation with 
that of DMO. The mass function is suppressed in the hydro- 
dynamic run by ~ 10% at 10 12 /i _ 1 M q , increasing to 20% at 
10 9 H~ 1 Mq. This agrees qualitatively with the findings from 
the OWLS simulation (Velliscig et al. 2014). The two-point 
correlation function for the dark matter particles in the two 
simulations is similar on large scales. On small scales, the 
two-point correlations are larger for the hydrodynamic run, 
corresponding to an increased dark matter power spectrum 
at large k values. 

Identifying matching subhalos in the two simulations 
enables us to directly compare the distribution of axis ratios 
and shape alignments to understand the effects of baryonic 
physics. The fraction of matched subhalos decreases as we go 
to lower masses, but is above 90 per cent for the mass range 
used in this work. These shapes are also used to calculate 
and compare the intrinsic alignment two-point statistics. 

We measured the shapes of dark matter subhalos in 
MBII and DMO simulations as a function of radius, by 
only using particles within a certain distance from the cen¬ 
ter based on a fixed fraction of the subhalo ifeoo- We an¬ 
alyzed the distributions of axis ratios and alignments of 
shapes measured within 0.2i?2oo, O. 6 R 200 , and I.OR 200 . In 
both simulations, we found that the axis ratios increase with 
the distance from center, which agrees with previous re¬ 
sults in iV-body simulations (Schneider, Frenk & Cole 2012). 


We also found that the dark matter subhalo axis ratios in 
the MBII simulation are higher (rounder) than those in the 
DMO simulation at all mass ranges due to the effects of 
baryonic physics, in agreement with previous findings (e.g., 
Kazantzidis et al. 2004; Abadi et al. 2010; Deason et al. 
2011). The galaxy stellar components in MBII have smaller 
axis ratios than the dark matter subhalos in both simu¬ 
lations, with the fractional difference being larger for the 
minor-to-major axis ratio, s. 

The degree of alignment of the stellar component in 
MBII with the dark matter component in DMO is larger in 
subhalos of high mass and decreases at lower masses. This 
trend is qualitatively similar to that of the alignment of the 
stellar component with the dark matter component in MBII. 
In subhalos of higher mass, the shapes of the dark matter 
component in MBII and DMO simulation are well aligned 
with each other. This alignment is stronger than the align¬ 
ment of the dark matter shape in MBII with that of stellar 
component in the same simulation. Comparing the misalign¬ 
ment of the stellar component with shapes of dark mat¬ 
ter component measured within different radial distances in 
MBII, we find that the misalignment angles increase as we go 
to larger radii within the dark matter subhalo. However, we 
do not notice a significant change in the misalignment of the 
stellar component with the shape of dark matter component 
in DMO measured for different radii. Based on our results 
from axis ratios and misalignment angles, we concluded that 
when mapping galaxy alignments from hydrodynamic sim¬ 
ulations onto subhalos in dark matter-only simulations, it is 
not useful to measure the shape of dark matter component 
at smaller radii in order to trace the shape and orientation 
of the stellar component. Hence, for the comparison of two- 
point statistics, we only consider the shapes of dark matter 
subhalos obtained using all the particles in the subhalo. 

Using the shapes of matched subhalos, we compared the 
intrinsic alignments two-point statistics (ellipticity-direction 
correlation, or ED, and the projected-shape density corre¬ 
lation, w§+) of the dark matter and stellar matter in MBII 
with that of the dark matter component in DMO simula¬ 
tion. The ED correlations of the dark matter component in 
MBII and DMO agree on large scales. On small scales, the 
ED correlation function decreases in MBII due to a change 
in the dark matter profile caused by baryonic physics. For 
the stellar component, the ED correlation is smaller on all 
scales due to the misalignment of the stellar component with 
the dark matter component. This corresponds to a frac¬ 
tional difference ranging from ~ 30 — 70% for scales around 
0.5—5/i _1 Mpc, and the decrease is larger in subhalos of lower 
mass. The w$+ correlation function for the shapes of dark 
matter component in MBII are smaller when compared to 
that of the DMO simulation due to the smaller values of el- 
lipticities, since the shapes are rounder in MBII. For the stel¬ 
lar component, we find that the w$+ is comparable on large 
scales and small scales with the W 5 + of the dark matter com¬ 
ponent in DMO. However, for scales around 0.5 — 5/i _1 Mpc, 
w§+ is still smaller for the stellar component by ~ 30 — 40%. 
At intermediate scales, we find a transition from the 1-halo 
to 2-halo term that causes a decrease in w$+ computed us¬ 
ing stellar shapes, but this feature is not clearly evident for 
the w§+ of the dark matter shape in MBII. 

Our results in this paper suggest that the scatter in the 
distribution of axis ratios of the dark matter subhalo shapes 
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Figure 14. Histograms of misalignment angles of 3D shapes of dark matter subhalos in the DMO and MBII simulations (with respect 
to the shape of the galaxy in MBII) when measuring the subhalo shapes within different radii (O. 2 R 200 , O. 6 R 200 , and 1 .OR 200 ) in the 
mass bins Ml (left), M2 (middle), and M3 (right) at z = 0.06. 



r{h l Mpc) r [h 1 Mpc) rf/T 1 Mpc) 


Figure 15. Comparison of the ED correlation function, cjg(r), for the dark matter subhalo and stellar matter components in MBII with 
respect to that for dark matter subhalos in the DMO simulation, computed separately for Ml (left), M 2 (midddle), and M3 (right). 



Figure 16. Comparison of the projected density-shape correlation, u>s+ (r p ), for the dark matter subhalo and stellar matter components 
in MBII with respect to that for dark matter subhalos in the DMO simulation, computed separately for Ml (left), M2 (midddle), and 
M3 (right). 
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Figure IT. Comparison of w§+(r p ) (projected shape-density correlation function without ellipticity weighting), for the dark matter 
subhalo and stellar matter components in MBII with respect to that for dark matter subhalos in the DMO simulation, computed 
separately for Ml (left), M2 (midddle), and M3 (right). 


in the DMO simulation is large compared to that of the stel¬ 
lar component in MBII, with significant misalignment in the 
orientation of shapes. However, the alignments between the 
galaxies in MBII and the corresponding matched subhalos in 
the DMO simulation are still significant compared to a uni¬ 
form random distribution. In future work, we will use these 
measurements to map the intrinsic alignments of the stellar 
matter component in hydrodynamic simulations onto dark- 
matter-only simulations. This is an important step in pro¬ 
ducing iV-body based mock catalogs that have realistically- 
complicated intrinsic alignments for tests of weak lensing 
analysis methods in future surveys. 
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